The Cronin Effect, Quantum Evolution and the Color Glass Condensate 
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We show that the numerical solution of the classical SU{3) Yang-Mills equations of motion in 
the McLerran-Venugopalan model for gluon production in central heavy ion collisions leads to a 
suppression at low pt and an enhancement at the intermediate pt region as compared to peripheral 
heavy ion and pp collisions at the same energy. Our results are compared to previous, Color 
Glass Condensate inspired calculations of gluon production in heavy ion collisions. We revisit the 
predictions of the Color Glass Condensate model for pA (dA) collisions in Leading Order and show 
that quantum evolution-in particular the phenomenon of geometric scaling and change of anomalous 
dimensions-preserves the Cronin enhancement of pA cross section (when normalized to the leading 
twist term) in the Leading Order approximation even though the pt spectrum can change. We 
comment on the case when gluon radiation is included. 

PACS numbers: xxx 



I. INTRODUCTION 

The Relativistic Heavy Ion Collider (RHIC) has 
opened a new frontier in high energy heavy ion collisions. 
The first data from RHIC which showed a large suppres- 
sion in the ratio of produced hadrons in AA and pp col- 
lisions Ql has created much excitement in the heavy ion 
community. It has led to intense theoretical and experi- 
mental work in order to understand and characterize the 
outcome. More recently, hadron spectra in dA collisions 
at RHIC were measured Q in order to clarify the role 
of and distinguish between initial state and final state 
(plasma) effects. Whether the observed suppression of 
hadronic spectra Q and disappearance of back to back 
jets in central Au-|-Au collisions, as well the large el- 
liptic flow [3, in heavy ion collisions at RHIC is due 
to the Quark Gluon Plasma is still in need of further 
experimental investigation. Nevertheless, the recent re- 
sults from the dA collisions at RHIC and the apparent 
lack of strong initial state effects in mid-rapidity and at 
high Pt appear to necessitate the presence of final state 
interactions in the partonic matter created in heavy ion 
collisions at RHIC @. 

Even though the Cronin effect (the observation, at 
fixed target experiments, that the ratio of pA and pp 
cross sections, scaled by the number of collisions, is above 
unity at some intermediate pt while below one at low 
Pt) is likely small in high energy heavy ion collisions as 
compared to parton energy loss effects, it is one of the two 
main nuclear effects, along with shadowing, expected in 
high energy pA (dA) collisions. Since pA collisions were 
proposed in order to clarify the role of initial and final 
state (medium) effects in high energy heavy ion collisions, 
it is extremely important to have a firm understanding of 
the physics of shadowing and the Cronin effect in order to 
have a precise understanding of the role of parton energy 
loss effects in heavy ion collisions. 



The Cronin effect in high energy pA and AA collisions 
has been the subject of renewed theoretical interest re- 
cently In this brief letter, we show that the nu- 
merical solution of the classical Yang-Mills equations of 
motion [TlL [T3. [T3 | in the McLerran-Venugopalan model 
[TH does indeed include the Cronin effect. We point 
out the differences between this numerical approach and 
other saturation inspired models [T^lisll which led to the 
absence of the Cronin effect in these models. We re- visit 
some of our earlier results for pA collisions and show that 
quantum evolution and the so called geometric scaling 
phenomenon |l9l l20l | (plus change of anomalous dimen- 
sion) preserves the Cronin enhancement even though its 
magnitude and peak may change. 



II. THE CRONIN EFFECT IN AA COLLISIONS 

In the McLerran-Venugopalan model, the classical 
Yang-Mills equations of motion arc solved in the pres- 
ence of external sources of color charge. These color 
charges can be thought of as the high x quarks and glu- 
ons in the wavefunction of a nucleus and are Gaussian 
distributed with a characteristic scale A^. (In practice, 
As « Qs, the saturation scale.) In principle, A^ can 
be determined from nuclear gluon distributions. It is 
an external parameter in the calculations described here. 
Physical quantities are computed by averaging over the 
Gaussian distributed color charges. The details of these 
computations for the real time gluodynamics of nuclear 
collisions can be found in Refs. pdl Il4l| . 

Briefly, the numerical lattice calculations in Ref. 
impose color neutrality condition at the nucleon level and 
realistic nuclear density profiles. In our computations, 
we first solve the classical Yang-Mills equations on the 
lattice for the two nuclei before the collision. In the radi- 



ation gauge {x~^A 



'A~^ = 0), the initial conditions 
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of gauge fields A'^ for nucleus-nucleus collisions at r = 
can be obtained by matching the solutions in the space- 
like and time-like regions. Requiring that the gauge fields 



be regular at r 



0, D^,F^'' 



and D„+F^'+ = J+ for 



gives the boundary conditions at t = 0: 



A^{Q,xt) 



yll(0,a;T) + A*2(0,2;t), {i^x.y), 



±x^ 



[^U0,a;T),A*2(0,XT)] 



Here, A\ 2 are the pure gauge fields for two incoming nu- 
clei. Using these initial conditions, classical Yang-Mills 
equations are solved, assuming boost invariance, on a 2- 
dimensional lattice. The definition of the number distri- 
bution in the non-perturbative region is discussed in de- 
tail in In the transverse Coulomb gauge V_l • ^ = 0, 
it is 



N{k) = V(l0(fc)^l)(k(fc)l^), 

where 4){k) and 7r(fc) correspond to the potential and ki- 
netic terms in the Hamiltonian respectively. 

In Fig. n we plot i?cp, the ratio of produced gluons 
in head-on {h — Q fm) and in peripheral (6 = 11 fm) 
Au-Au collisions normalized by the number of collisions 
-^coii, for an SU(3) gauge theory. A'^coii is computed self- 
consistently in our framework and agrees, for instance, 
with Ref. |1S]. The ratio Rcp is below unity at low 
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FIG. 1: Rcp from the McLerran Venugopalan model for an 
SU(3) gauge theory. Here A^o ~ 2 GeV, where A^q is the 
color charge squared per unit transverse area in the center 
of each nucleus. (The value of As averaged over the entire 
nucleus is smaller ~ 1.4 GeV.) This result is obtained for a 
256 X 256 lattice. 

pt and above unity at intermediate pt- This shows that 
the original McLerran- Venugopalan model indeed has the 
Cronin enhancement. One can in principle show that 
this ratio goes to unity at high pt, but this is numerically 
intensive and not shown here. 

Instead, for simplicity, we show in Fig. [21 results from 
a computation of the real time evolution of an SU(2) 



MV 

-color neutral I 
-color neutral II 



gauge theory. There is no qualitative difference between 
the SU(2) case and the SU(3) case. Our results for the 
SU(2) gauge theory are obtained for larger (512x512) 
lattices relative to the (256x256) lattices in the SU(3) 
case. Plotted in Fig. |21 is the ratio of the pt distri- 
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FIG. 2: Raa for an SU(2) gauge theory. Raa is the ratio 
of the Pt distribution of gluons for Au-Au collisions (Aso = 2 
GeV) divided by p-p collisions (Aso = 0.2 GeV) and normal- 
ized by the ratio of their asymptotic values at large pt- Here 
MV denotes the McLerran- Venugopalan model with color 
neutrality imposed globally; Color Neutral I & II impose color 
neutrality on each configuration at the nucleon level-see text 
for discussion. 

bution of gluons from Au-Au collisions (with A^o = 2 
GeV) to that in p-p collisions, normalized to the ratio of 
their asymptotic values (TVcoii) at high pt- The p-p re- 
sults here are taken to be the lattice results for the very 
small value of A^o = 0.2 GeV. The lattice result in the 
latter case is equivalent to the perturbative tree level re- 
sult up to very small pts. The reason Raa deviates from 
unity is due to the multiple scattering ("Cronin") effect 
illustrated in Fig. lathis point will be discussed further 
shortly. The three curves in Fig. [21 correspond to the fol- 
lowing: MV denotes the McLerran- Venugopalan model 
with color neutrality imposed only globally over the en- 
tire nucleus in each configuration. Color Neutral I & II 
correspond to the more stringent condition where color 
neutrality is imposed on the scale of a nucleon. The for- 
mer corresponds to the case where monopole component 
of the nucleon color charge is subtracted while in the 
latter case, both monopole and dipole components are 
subtracted. As has been noted previously in Ref. [T^ . 
the effect of these more stringent color neutrality condi- 
tions is to induce a power law dependence at low pt for 
the correlator of color charges in momentum space. This 
dependence, albeit non-perturbative in origin, is similar 
to that induced by perturbative color neutrality which 
arises from the color screening of saturated gluons gener- 
ated in the quantum evolution [T9l | . The softening of the 
Cronin enhancement can therefore be understood as re- 
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suiting from color screening or "shadowing" of the initial 
giuon distributions. 




FIG. 3: Schematic diagram of mono-jet giuon production am- 
plitude in kt factorized form (left figure) vs the non-factorized 
contributions arising from solutions of classical Yang-Mills 
equations (right figure). 

In a nice paper, Kovchegov and Mueller have shown 
that the choice of gauge can determine whether interac- 
tions are "initial state" or "final state" effects . In our 
numerical computations (carried out in = gauge) 
initial state as well as some final state interactions are 
included. We qualify the latter because due to the strong 
expansion of the system with time, giuon field strengths 
decrease very quickly. High pt modes (pt > Ago) in par- 
ticular arc linearized at very early times and therefore, 
suffer no further final state interactions. On the other 
hand, soft modes still interact with each other for a long 
time as shown in Ref. Therefore, the lattice result 

of Rc p > 1 in Fig. ^ can be regarded as the Cronin 
enhancement resulting from these re-scatterings. The in- 
teraction between hard and soft modes for occupation 
numbers / < 1 remains an outstanding problem. Naively, 
these interactions are of higher order in the classical ap- 
proach but at late times, with the classical approach 
breaking down, they become competitive. Despite recent 
progress |22| . a full treatment of this problem is lacking. 

We now compare our treatment to other approaches in 
saturation based models. Since there is no known com- 
plete analytical solution for the collision of two nuclei 
psL 01 J a simplified model was used in . In this "fct 
factorized" approach, the cross section for giuon produc- 
tions is given by 

^da AttNc 1 r* „2 , / ,2^ , / . , ^2^ 

d-^p -IPt J 

where (j)A{x^pt) is the unintegrated giuon distribution 
function of a nucleus 

GAix,Q^) = J dk^Mx^h). 

This equation describes mono-jet production in the kt 
factorized form and can be visualized as in Fig. 01 Ra- 
diation of additional jets are down by powers of a^. The 




FIG. 4; Mono-jet production cross section in kt factorized 
form. 

connection to giuon saturation and the Color Glass Con- 
densate is phenomenological and comes from the assumed 
form of the unintegrated giuon distribution <j)A{x,kt) in 
the transverse momentum region ^ Q^. In the per- 
turbativc region where fcf ^ 

cj)A{xukl)^l/k^t (1) 

while in the saturation region where k^ <^ Qs 

^A{xi,kl) ^l/as. (2) 

In these two asymptotic forms of the unintegrated 
giuon distribution are then matched at Qs- In this ap- 
proach, there is only a single scattering in the kt > Qs 
region and multiple scatterings at scales just above Qs 
are not included. Since the non-perturbative input in the 
form of the "giuon liberation factor" is included from the 
numerical lattice simulations this approach likely 

includes phenomenologically the physics of the CGC for 
global observables such as the centrality and rapidity de- 
pendence of observables (further discussed in the next 
paragraph) . 

The same fct-factorized formalism is considered in 
Ref. . The authors however additionally consider the 
effects, at moderate pt's, due to the change in the unin- 
tegrated giuon distributions arising from quantum evolu- 
tion. In particular, they take into account the change in 
the distributions due to the modification of the anoma- 
lous dimensions in the evolution. This change leads to 
a scaling with A^part as opposed to A^coU scaling. (Note 
that these modifications do not affect the prior results 
of Ref. ^3 since the multiplicities, at the pt's at which 
they occur, are rather small.) However, the Cronin en- 
hancement for Pt distributions is missed in 0, 0| . On 
the other hand, in our numerical computations, all re- 
scatterings are included at the classical level leading to 
the Cronin enhancement. Note though that quantum 
evolution is absent in our lattice calculations-its effects 



4 



are only implicitly included through the magnitude of the 
saturation scale. 

In addition to their importance for the Cronin ef- 
fect, the non-factorizcd contributions illustrated in Fig.|2| 
likely play a significant role in determining the total gluon 
multiplicity as well. This can be understood by compar- 
ing the numerical results for the gluon liberation coeffi- 
cient /at 0,^3 which give = 0.3 — 0.5, in contrast 
to analytical results for this coefficient that include only 
the factorized contributions and is roughly 3 to 4 times 
larger Q. 

A possible way of including this enhancement in the kt 
factorized form through quantum evolution is discussed 
in psf . A similar approximation of 0a (a^, kt) in the con- 
text of pA collisions also leads to a lack of Cronin 
enhancement in the cross section. It would be interest- 
ing to include quantum evolution effects in addition to 
the Cronin effect in our numerical work to quantify the 
importance of the energy loss contribution. This work is 
in progress and will be reported shortly |2^. Below, we 
consider pA scattering and explicitly show that inclusion 
of higher order scatterings leads to the Cronin effect. 



III. THE CRONIN EFFECT IN pA COLLISIONS 

Even though the Cronin effect in a high energy heavy 
ion collision is quite small at high pt as compared to 
the parton energy loss effects, it is quite important in 
a proton (deuteron) nucleus collision since energy loss 
effects in the hot medium produced in a heavy ion col- 
lision, are absent in a pA collision. In this section, we 
re-visit the Color Glass Condensate results for a quark 
or gluon scattering on a nucleus Q and show that the 
recent phenomenon of geometric scaling and changing of 
anomalous dimension does not affect our conclusions as 
long as we divide our pA cross section by its leading twist 
limit (rather than pp) . 

Scattering of a quark or gluon on a nucleus using the 
Color Glass Condensate formalism and its relation to DIS 
was discussed in |^ . It was shown that the scattering 
cross section is given by 

<Pptdp-<Pbt " (2^)2 ^^"^ P ' 

X J dVte^P* '-*[2-(7dipoic(a;,rt,60](3) 

with 

adipoic(a;,r„ 50 ^ ^Tr (l - Vih + ^-^)V\bt - ^))^(4) 

where q" [p~ ) is the longitudinal momentum of the in- 
coming (outgoing) quark, with a similar equation for 
gluon scattering. This is the multiple scattering general- 
ization of quark gluon scattering in pQCD and unlike the 
leading twist (single scattering) result, is finite as P( — > 



due to higher twist effects. In addition, defining 

^{x.ptM)^ y'd^r.e^P'-'^' [2-adipoie(x,rt,6t)] (5) 

and using the fact that Cdipoic (a;, rt = 0,6*) = 0, we see 
that 7 satisfies the following sum rule 

J d^pa{x,pt,bt)^ 2 {27r)\ (6) 

for fixed bt- It is clear from this sum rule that if the 
cross section in © is suppressed at low pt, it must be 
enhanced at high pt in order to compensate for the low pt 
suppression so that the sum rule holds |30l |. The Cronin 
peak moves to higher pt as one goes to higher energy (or 
more forward rapidities as well as more central collisions). 

The effects of quantum evolution on our classical 
results can be investigated using the Rcnormalization 
Group equations 1271 . and in particular the Balitsky- 
Kovchegov (BK) [23 equation for the dipole cross sec- 
tion in the large Nc limit. It is straightforward to apply 
the BK equation to the qA (or gA) scattering cross sec- 
tion in order to prove that quantum evolution preserves 
the sum rule in ©. For our purposes, it suffices to no- 
tice that adipoioix, rt = 0, 6t) = at any x (energy) so 
that the sum rule is preserved by quantum evolution and 
therefore, the Cronin enhancement is still present. Nev- 
ertheless, the Pt spectra will look different at different 
energies since the location of the Cronin peak as well as 
its magnitude changes as one goes to higher energies (see 
for example, Figure 4 in the last paper of [^.). In ad- 
dition, the change of the anomalous dimension 7 could 
modify the magnitude of the Cronin effect. Nevertheless, 
since the BK equation preserves the sum rule in @ at 
the partonic level, this modification in the high pt region 
will be compensated by an analogous modification in the 
low Pt region so that the sum rule is not violated. 

One should keep in mind that at RHIC and for mid 
rapidity and high pt processes (such as the suppression 
of hadron yields) the effective x of the partons is quite 
large. For example, for pt '^5 — 10 GcV in mid rapidity, 
the X range is ~ 0.05 — 0.1. The only potential evidence 
for gluon saturation and geometrical scaling in hadrons is 
from HERA where saturation models work very well only 
for X < 0.01 and fail at higher x. In the case of nuclei, 
this effective x may be slightly higher but can not be 
too much higher since otherwise, strong gluon shadowing 
would manifest itself in the F2 structure function which 
shows no strong shadowing effects at the x range 0.05 — 
0.1 and as a matter of fact, shows anti shadowing! The 
Color Glass Condensate model is an effective theory of 
QCD at small x only and likely breaks down at the high 
x's relevant for high pt processes in mid rapidity at RHIC. 
(This situation will improve as one goes to higher collision 
energies or stays in the low pt region.) 

Moving to more forward rapidities (in the projectile 
fragmentation region) will make applications of the Color 
Glass Condensate model more reliable since smaller val- 
ues of X in the target are probed. The saturation scale 
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is larger, rendering weak coupling methods more reliable. 
In addition, the contribution of high x region to the cross 
section will be less important than that in the mid ra- 
pidity region. Whether the most forward rapidities ac- 
cessible at RHIC are large enough for the physics of gluon 
saturation to be the dominant physics remains to be seen. 

IV. SUMMARY 

We considered in this note gluon production in heavy 
ion collisions using the numerical simulations of the Color 
Glass Condensate. We showed that the ratio of central 
to peripheral cross sections (or equivalently, the ratio of 
central to pp cross sections), normalized by the number 
of collisions, shows the Cronin enhancement at high pt 
as well the suppression at low pt- We discussed other 
gluon saturation and Color Glass Condensate inspired 
models and commented on the absence of the Cronin ef- 
fect in these models. We also considered pA collisions and 
showed that, at the partonic level and in leading Order 
in as, quantum evolution with energy and the change of 
anomalous dimension preserves the Cronin enhancement 
(when normalized to its leading twist term rather than 
pp) due to a sum rule satisfied by the dipole cross section 
even though the magnitude and location of Cronin peak 



is energy dependent. When divided by the pp cross sec- 
tion as done experimentally, our ratio will have suppres- 
sion at high Pt in agreement with the results of Ref . . 

After this work was completed, we were made aware 
of similar work by Kharzeev, Kovchegov and Tuchin |29l | . 
Although there is some overlap in our discussion of the 
pA case (their focus here being on gluon production), 
they do not explicitly consider the Cronin effect for AA 
collisions as we have. 
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